Stationary phase corrections in the process 
of bosonization of multi-quark interactions 



A. A. Osipovt 1 , B. Hillert 2 , J. Moreira t3 , A. H. Blin^ 4 

t Centro de Fisica Teorica, Departamento de Fisica da Universidade de Coimbra, 

3004-516 Coimbra, Portugal 



Abstract 

The functional integration over the auxiliary bosonic variables of cubic 
order related with the effective action of the Nambu - Jona-Lasinio model 
with 't Hooft term has recently been obtained in the form of a loop expan- 
sion. Even numbers of loops contribute to the action, while odd numbers of 
loops are assigned to the measure. We consider the two-loop corrections and 
analyse their effect on the low-lying pseudoscalar and scalar mass spectra, 
quark condensates and weak decay constants. The results are compared to 
the leading order calculations and other approaches. 



PACS number(s): 12.39.Fe, 11.30.Rd, 11.30.Qc 



On leave from Joint Institute for Nuclear Research, Laboratory of Nuclear Problems, 
141980 Dubna, Moscow Region, Russia. Email address: osipov@nusun.jinr.ru 

2 Email address: brigitte@teor.fis.uc.pt 

3 Email address : j moreira@ teor . fis . uc.pt 

4 Email address: alex@teor.fis.uc.pt 



1 



1. Introduction 



The large distance dynamics of QCD is dictated to a great extent by the 
spontaneous symmetry breaking of chiral symmetry QjJ The Nambu - 
Jona-Lasinio (NJL) model of fermionic fields j.3] suggests that the dynami- 
cal mechanism for such breaking be in analogy with the Ginsburg - Landau 
theory of superconductivity Numerous studies |3]-[7j have been per- 
formed since that time with the respective effective mesonic action derived 
from four-quark interactions of the NJL type. During these years the res- 
olution of Ua(1) problem has been found and, in particular, the relevance 
of the U(3)l x U(3)r chiral symmetric NJL model combined with the six- 
quark 't Hooft flavour determinantal interaction (NJLH) [S| for low-energy 
phenomenology of mesons was noted [HI- [121 ■ The explicit breaking of the 
unwanted axial symmetry by the 't Hooft determinant is motivated 

by the instanton approach to low-energy QCD [H], [T3j . 

Originally written in terms of fermionic degrees of freedom, the NJLH 
model has been widely explored at mean field level with Bethe-Salpeter and 
Hartree-Fock techniques applied to quark-antiquark scattering in its various 
channels of interaction j§] , [TH Ho] . 

In parallel, functional integral methods have been used to obtain the 
Lagrangian in bosonized form ^U] 5 jTB]- . The bosonization gives rise to a 
doubling of the mesonic auxiliary fields, of which one set has to be integrated 
out. This latter, in the presence of the 't Hooft interaction, involves a term 
of cubic order, which cannot be integrated out exactly. In [10] the leading 
order stationary phase approximation (SPA) was calculated. At this order 
the effective potentials obtained with both methods coincide |17j . 

Given its success in describing a large bulk of empirical data, the question 
arises whether corrections to the leading order SPA result are small. By em- 
barking in this task, a series of startling facts came across our investigations 
BB|: 

(i) The stationary phase equations which one obtains in the NJLH model 
have more than one root (critical point). Only one has a regular behaviour 
in the limit k — > of the six-quark coupling, the others are singular. The 
rigorous SPA treatment requires taking into account all critical points, which 
give rise to an unstable vacuum for the theory. 

(ii) The result obtained in jTUj corresponds to the regular root contri- 
bution. It is an approximation which leads to an effective potential with a 
well separated local minimum, which approaches smoothly the stable NJL 
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vacuum as k — > 0. Such a local minimum is probably a good ground for phe- 
nomenological estimates, at least all known calculations made in the NJLH 
model are based on this approximation. It has been shown recently fH] 
that eight-quark interactions stabilize this vacuum state, opening the way to 
justify this approach theoretically 

(hi) Two expansions of the effective action have been considered: the 
perturbative series in k and the loop expansion. Both of them were never 
studied beyond the leading order. It is tacitly assumed that next to the 
leading order corrections are small, although this fact has never been proven. 

In this paper we quantify the two-loop order contributions to the La- 
grangian derived previously by studying their impact on the mass spectrum 
of low-lying mesons. We show that the effect is of the order of a few percent 
compared to the leading order masses, improving them slightly. We think it 
is rather safe to conclude that the loop expansion is rapidly converging, at 
least for the mass spectra. 

The paper is structured as follows. In the next section we collect the 
essential information needed to extract the linear and quadratic terms which 
contribute to the gap equations and mass terms respectively. In section 3 we 
write out the expressions for the gap equations, masses, weak decay constants 
and condensates. Section 4 contains the numerical results and discussion. 
Conclusions are presented in section 5. 

2. The "tandem" Lagrangian 

To be self-contained and define the notation we review the main ingredi- 
ents of our model calculations. The underlying multi-quark Lagrangian is 
bosonized in a two step (tandem) process, in which a semi-bosonized func- 
tional, quadratic in the fermionic fields, and another functional, depending 
only on the auxiliary bosonic variables, can be dealt with separately. The 
integration over the quadratic fermionic degrees of freedom is formally exact 
and is calculated using a generalized heat kernel method. We start by pre- 
senting how these two types of functionals emerge and how we obtain and 
calculate the loop expansion we are after. 
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2. 1 The stationary phase contribution 

We consider the fermionic Lagrangian 

Cnjlh = qiiYd^- m)q + £ NJL + £ H , (1) 

which contains the NJL four-quark vertices of the scalar and pseudoscalar 
types 

G 

£n.jl = — [{qKq) 2 + (qh^ a q) 2 \ (2) 

and the six-quark 't Hooft interaction |5] 

C H = «(det qP L q + det qPaq), (3) 

where m is the diagonal current quark matrix for quark fields with Nf = 3 
flavours and N c = 3 colours. In eq.(J2J) A a , a = 0,1. ..8, are the normalized 
(tr\ a \b = 25 a b) matrices in flavour space. The explicit form of these U(3) 



hermitian generators is A = y2/3 and A a for a are the usual Gell- 
Mann matrices. The positive coupling G, [G] = GeV -2 has order G ~ 1/N C . 
In (jHJ) the negative coupling k of dimension [k] = GeV -5 has the large N c 
asymptotics k ~ 1/N C f ■ Therefore Cnjl dominates over Ch at large A^" c . 
The matrices Pl,r — (1 T 7s)/2 are projectors on the chiral states and the 
determinant is over flavour indices. 

We are assuming that the quark fields transform like the fundamental 
representations of the global U(3)l x U(3)r chiral group, i.e., 

Sq = i( a + 7B/9)g, 5q = -iq(a - 75^), (4) 

where the parameters of the infinitesimal transformations are chosen as a = 
®a^a, P = Pa^a- One now observes that 

5Cnjlh = iq[[a, m] - ^{P, m}) q + 2iV6p n (det qP R q - det qP L q) . (5) 

The global chiral symmetry is broken explicitly by the current quark mass 
term and the U(1)a axial symmetry is broken too due to the 't Hooft inter- 
action. 

The functional integral in bosonized form is derived in jlOj . and has the 
form 



Z = J Y[VUA'DqVq exp J d 4 xC q (q,q,a,(, 

+00 

x J Y[VR A exp^ J d A xC r (U, A;R)j , (6) 
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where 



C q = qih^ - M — a — * 75 0)g, (7) 
C r = R A (U A + A A ) + ^ R\ + | $ ABCj R Aj R Bj R c . (8) 

with a cubic polynomial in the fields R A in the exponent. The notation is 
as follows [IE]: R A = (R a ,Ra) = (s a ,Pa) and U A = (n^ILj = (cr a ,<f>a) are a 
very compact way to represent two sets of auxiliary bosonic variables, each 
containing a scalar s a (cr a ) and a pseudoscalar p a (0 a ) nonet. The indices 
(a, d) run from to 8 independently in flavour space. We also define the 
related quantity A A = (A a , 0) = (M a — m a , 0). 

The external scalar fields o = cr a A a have been shifted a a — > a a + M a by 
the constituent quark mass M a , so that the expectation value of the shifted 
fields in the vacuum corresponding to dynamically broken chiral symmetry 
vanish. The vacuum expectation value of the "unshifted" scalar field 

(a) = M a X a = diag(M u , M d , M.) (9) 

gives the point where the effective potential of the model V((cr)) achieves its 
local minimum. The corresponding condition is known as the "gap" equation. 
It eliminates tadpole graphs and determines the values of constituent quark 
masses as functions of the model parameters and of the cutoff A. The case 
m u ^ md 7^ m s corresponds to the most general breakdown of the SU(3) 
flavour symmetry, giving M u 7^ Md 7^ M s . In this way the ground state 
of the system includes effects of the explicit symmetry breaking. We will 
assume in the following that m u = rrid. 

The variables a a and (j) a must be replaced by the physical scalar and pseu- 
doscalar states a% h and <^ determined through the appropriate normaliza- 
tion of their kinetic terms. Note that these terms, as well as other important 
contributions to the meson masses and interactions of the effective mesonic 
Lagrangian, are obtained as a result of integration over the quark fields in 
Z. One has 

<r a , <\> a = gaf, gctf a h (10) 

at leading order of the heat kernel expansion of the effective mesonic action 
(see next section). The quark-meson coupling g, being a function of param- 
eters of the model, fulfills in addition the Goldberger - Treiman relation at 
the quark level: g = M u /f n . Combining this relation with (JSJ) and (jl(J|) one 
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finds the well-known linear sigma model result (201 



K h ) = U- (ii) 

Finally, the three index coefficients &abc are defined as 

3 

®abc = -® ab6 = —A abc , $ o66 = $ &b6 = 0, (12) 

obeying 

®abc6bc = 0. (13) 

The totally symmetric constants A abc are related to the flavour determinant, 
and equal to 

A a bc — 7y Cijk e rnnl{^a)im{^b)jn{^c)kl- (14) 

Now the functional integral over the auxiliary variables R A in (jUJ) 

+00 

Z[U, A] = J ]JVR A exp (i J d 4 x£ r (n, A;R)) , (15) 

— oo A 

can be written in the form |TS] 

Z[U, A] ~ exp (i J d 4 x£ s t 
+00 

x J ]JVR A exp (*-J A^x£! A qRaRb 

-00 A 

x E - ^ $ abc / ^xR A R B Rc) . (16) 

n=0 n - \ 3. ■/ / 

Here £ st is the stationary value of the Lagrangian £ r (U,A;R) associated 
with the regular critical point, around which the effective Lagrangian has 
been expanded. The barred fields indicate that they are shifted with respect 
to their stationary values R = R — R st . We denote by C" AB the coefficient of 
the expansion to second order in the fields. The expansion stops exactly at 
the third order C" ABC = k&abc/3\, whose exponent is represented here as an 
infinite series. 

The terms with odd values of n do not contribute in (fTBj) . because the 
corresponding functional integrals over R A equal to zero. The first term, 
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n = 0, sums all tree diagrams of a perturbative series in powers of the 
coupling k resulting in £ st . The n = 2 term represents the first non- leading 
correction to the effective Lagrangian (for details we refer to where we 
show that this correction can be associated with a "two-loop" contribution 
of quantum auxiliary bosonic fields) 



£ 



eff 



£*.+ ^ 



3k 2 M 



2tt J 32N(N + 2)(N + 4) 



(17) 



The stationary Lagrangian reads to cubic order in the fields ITT] 



-St 



h a cr a + -h^a a a b 



O (field 



abc 



2 

, 

'abc 



1 hab4 > a4>b 



h {2) +h (3) ' 



>bVc 



(18) 



The two-loop corrections, contained in the second term of (|17|h consist of 
a rather intricate dependence 



M 



tr £"- 1 ) 3 + 6tr tr(£"- 1 ) 2 + 8tr (C'^j 



(19) 



on flavour traces of powers of C" , which is the inverse of the real and 
symmetric N x N matrix (N = 18) 



St I 



( 



3/Sj 3/x 
GS a b + — A a b c sl t ——A a b c p c st 
o 16 16 

_ 1~6 fcc ^* st * jr " ab ~ Yq abcS st 



(20) 



calculated at the stationary points s£. and p^.. These are expressed in in- 
creasing powers of the external fields a a , (fi a 



3 St 



+ habcdWcVd + h { 2 cd O h 4>c4>d + • • • 



hab<Pb + KL^bOc + h^ b ' cd a b cr c (j) d 
oiH-.mMd + ■■■ 



,(3). 



,(3) 



(21) 



(22) 



with h a % l,„ depending on A a and coupling constants (see Appendix). In par- 
ticular the coefficients have nonvanishing components for a = (0, 3, 8) and are 
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obtained, with h = h a X a = diag(/t u , h d , h s ), in the case of isotopic symmetry 
{K = h d ) as [17J 



AC 

Gh u + A u + — h u h s = 0, 
lb 



Gh s + A s + ^-h 2 u = 0. 
lb 



(23) 



In eq. ()17|) A denotes an ultra-violet cutoff associated with the stationary 
phase corrections to the functional integral over auxiliary bosonic fields. It 
is a free parameter to be fixed by phenomenology. 

To handle the new contribution Ai, we expand C" AB (R st ), which we ab- 
breviate from now on as £", to second order in the external fields a a , (fi a , 

£" = L + L 1 + L 2 + O(field 3 ), (24) 

and its inverse 

C"- 1 = L + L 1 + L 2 + O(Re\d 3 ). (25) 

This is all one needs to extract the relevant terms to the masses arising from 
the two-loop correction term. Here Li, (i = 0,1,2), denote the matrices 
which are constant, linear and quadratic in the fields respectively. The Li 
are constructed order by order, starting from the 0th order 

C'C"- 1 = L L = 1, (26) 

i.e. Lq = Lq 1 . 

The next terms are conditioned by this relation. Combining the first 
order Lagrangians and truncating at the linear fields 

(L + L 1 ){L + Lx) -> L L + L 1 L + LqIx = 1 (27) 

one gets the matrix Z 1; after using (J2HJ), 

L\ = —LqLxLq. (28) 

In a similar fashion one derives the matrix L 2 as 

L L 2 L + LqL 1 L 1 ) = — {LqL 2 Lq — LqLiLqLiLq) . (29) 
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Using Li in (|25|l and inserting in ()19|). one obtains 



trLs 



M 2 = 3 MtrLoJ (trLiJ + (trL 

+ 6 {tr (I ) [tr (ifj + 2 tr (l L 2 
+ 2 tr (Li) tr (l^) + (trL 2 ) tr (if}}, 



8 tr (Z Z? + L 2 Ll 



(30) 



where .M2 stands for the part of M. which contains only the second order 
terms in the fields a a ,(fi a . 

This expression is used in "Mathematica" [21] • Although the results, after 
evaluating of traces, are analytical they are very lengthy and not illuminating, 
and will not be presented here. However some structures are relevant for the 
low-energy theorems and the results will be encrypted in them, as shown in 
the section where we present the mass formulae. 



2.2 The heat kernel contribution 



It still remains to evaluate the functional integral over the quark degrees of 
freedom in eq.©. The Lagrangian C q is invariant under the chiral transfor- 
mations (@J) and the transformations 

5a = i[a,a + M] + 5<j) = i[a, 0] - {(3, a + M}, (31) 

induced by them for the external fields. All symmetry breaking terms have 
been absorbed in C r . This fact is of importance, since one can use then the 
generalized asymptotic expansion of the quark determinant |22[ 12*3*] . This 
method preserves the above mentioned symmetry at any order, taking into 
account the effects of the flavour symmetry breaking contained in the mass 
matrix M. Thus the corresponding part of the effective action can be written 

1 r 00 

In | detD\ = --— / d 4 x E ]T / i - 1 tr(6 i ), (32) 
327T J i=0 

where D = i^^d^ — M — a — is the Dirac operator present in C q , eq.(J7J), 
and the hi are generalized Seeley-DeWitt coefficients |22j . of which we show 
the first four for the case of SU{2)i x U(l)y flavour symmetry 

6o = l, 
&! = -Y, 



9 



Y 2 A 



It s 



Y 3 A 2 . A,, a . ... 1 



2 " 



h = -— + ^x 8 Y — ^ A 8 r 2 — (dr) 2 . (33) 

3! 6^3 2^3 12 v ; K J 

In the present case the background dependent structure Y is given by 

Y = i7„( V + ^5<9 M 0) + a 2 + {M, a} + 2 + z 75 [<7 + M, 0]. (34) 

We use the definition Ay- = M 2 — M 2 . In eq. (J23) the trace is to be taken 
over colour, flavour and Dirac 4-spinor indices and the regulator-dependent 
integrals I, are the weighted sums 

Ii = ~ (2J,(M 2 ) + J,(M 2 )) (35) 

with 

oo 

J,(M 2 ) = / ^-p(tA 2 ) exp(-tM 2 ). (36) 
o 

They are regularized with the Pauli-Villars regularization scheme [21] with 
two subtractions and one ultra-violet cutoff A 

p(tA 2 ) = 1 - (1 + tA 2 )exp(-iA 2 ). (37) 

We obtain, for instance, j2S] 

MM 2 ) — A 2 - M 2 In fl + ^) . ( 38 ) 

Jl(M »)=h,(l + ^)- F ^ S j. (39) 

Both of them are divergent in the limiting case A — > oo. Note that A does 
not need to be the same cutoff as A of eq.lfTTj). In the following we restrict 
our study to the two nontrivial terms, b\ and 62, i n the asymptotic expansion 
of In I det.D|. In this case only I and I\ are involved, related to the quark 
one-loop integrals of one- and two-point functions respectively, at zero four- 
momentum transfer. 
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3. Gap equations, condensates, meson spectra 



3.1 Gap equations and condensates 



The complete effective bosonized Lagrangian 



(40) 



comprises contributions from the heat kernel expansion to order 6 2; £hk, 
and from (|17|) . where C c stands for the two- loop corrections. We restrict to 
the case SU(2) I x U{l) Y symmetry, e.g. M u = M d ^ M s . The first two 
contributions remain the same as in the leading order calculations |12j . 

Equating the coefficient of <7j, % = (u,d,s) in eq. (}4*0|) to zero we obtain 
the gap equations 



h u + -^ M u [3/ + (M s 2 - M 2 u )h] + 2c u = 



h s + 



6tt 2 
Nc 
6tt 2 



M, 



3/ - 2{M 2 S - Ml)h 



+ 2c s 



0. 



(41) 



where c u and c s denote the corrections arising from C c . They depend on 
h u , h s , A, k. These equations must be solved self-consistently and in conjunc- 
tion with the stationary phase conditions The solutions Mi of (jUJ 
allow to calculate the condensates {uu) and (ss) (see eq.(|H4^)) 



4tt 2 



MiJ (M?) - miJ {m] 



(42) 



where we have subracted the contribution from the trivial vacuum jO]. Al- 
though they are structurally identical to the condensates calculated at lead- 
ing order, they encode the information of the correction terms q implicitly 
through Mj. 



3.2 Meson masses 



The expressions for the leading order masses, i.e. with C c put to zero, will 
not be repeated here. They were obtained in |12j . The correction mass terms 
can just be added to the leading order terms in their "raw" form, that is, 
as they are directly extracted from Chk, depending on Io,Ix integrals. To 
check the low energy theorems, one can then use the new gap equations, with 
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the correction terms q included, to eliminate these integrals. For example 
for the pion, </>j (J = 1, 2, 3), one has 

^H K (ml) = -^(3h + A„I 1 )^, (43) 

With "Mathematica" we are able to identify 

where Ui = nhi/(lQG). This connection between the pion mass correction 
to the gap equation correction term c u is crucial to guarantee the Goldstone 
limit. Indeed we obtain, after eliminating J , I\ from (}4*3*)) with help of the 
gap equations, 



m 



2 



° 2 ft i o ° u 



° 2 9 2 m u 

m * = M.G(l +.,.) ■ (46) 



where is structurally identical with the leading order pion mass; g 2 = 
4tt 2 /(N c I 1 ) renormalizes the pion fields to the physical fields (see eq.fJTUJ)) 
We used also that 

K 1 



A u G(l + u s ) 

which is a simple consequence of the stationary phase conditions (|2 
In an analogous way we are able to get for the kaon mass 



(47) 



02 _ g 2 (m u + m s ) 

m «- G(M u + M s )(l+u u ) ' ^ 

o 2 

where m K has the form of the leading order kaon mass. We used again (|23j) 
to obtain 

_ h u + h s = 1 

A u + A s G(l+u u ) ■ y } 
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Concerning the rj, 7/ corrections, we show here only some relevant prop- 
erties obtained with the help of "Mathematica" in the SU(3) limit 



(A<)c = (Am^) c = (Am 2 8 ) c , 
{Am 2 08 ) c = 0, 

(Am 2 ) c - (Am 2 88 ) c ± 0, (50) 

where, for instance, (Am 2 ) c is the contribution to the pion mass obtained 
from the Lagrangian C c . Therefore the correction to the flavour (0, 8) com- 
ponents follow the same patterns as in the leading order case [T2|, complying 
with the low-energy requirements: the octet member m 88 remains degener- 
ate with the pion and kaon, the mixing m^ 8 vanishes and in the chiral limit 
the correction to the mass of the singlet Wq is also non- vanishing, and will 
therefore contribute to the singlet-octet splitting. 

For the scalars we obtain in the SU(3) limit also that the corrections to 
the masses behave as 

(AM 2 J C = (AM 2 ) C = (AM 8 2 8 ) C , 



(AM 2 8 ) C = 0, 

(AM 2 ) C - (AM 8 2 8 ) C ^ 0. (51) 
3.3 Weak decay constants 

We use PCAC and the Gell-Mann-Oakes-Renner (GOR) j2i relation to ex- 
tract the condensates. From PCAC the weak decay constants are given as 

, M u M u + M s 

Jn = , JK = ~ • (52) 

9 2 9 

Using this and the mass relations for m w and m K , eqs.fHoT) and (jlHj) . one 
obtains the GOR equations with some model corrections of higher order in 
the current quark masses and from which one identifies the condensates 

m 2 J 2 = m u {h u + 2c u ) (l + ^) = -2m u (0\uu\0) (l + ^ 
-2 f 2 _ 1/ , w, , „ , 1. , f, , m u + m s 



m KfK = -^( m u + m s )(h u + 2c u + h s + 2c s ) (l + ^ ^ 

= -\{m u + m s )(0\uu + ss\0) (l + a' ) ' (53) 
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Finally, by using the gap equations (jlTj) in (|53|) and expressing I , ii 
through J (M 4 2 ), J\(Mf) with we obtain the condensates as 

(0|^|0) = -£| M^M 2 ) + 0(J 2 ), (54) 

where the 0(J2) terms are neglected, to conform with the truncation of the 
heat kernel series. We recall that the 0{J2) emerge from a property of the 
generalized heat kernel series in which differences of Jfc(M 2 ) — Jk(M^) are 
expressed as an infinite series involving Jk+i, I > [22] . 



4. Numerical results and discussion 

There are six parameters in the model, m u , m s , G, n, A, A. To see the effects 
of the new contribution, proportional to the cutoff A, we compare pairwise 
in the sets (a,b), (c,d) and (e,f) of Tables 1-3 the results calculated with 
A = and A / 0, keeping the remaining input unchanged. In this way, 
within each pair of sets, a running value of A between the indicated ones, 
will interpolate smoothly between the calculated observables shown. In a-d 
we fix four parameters through the pseudoscalar sector, m n , m K , f n , f K , and 
adjust k through the quark condensate (uu). In sets (e,f) five parameters are 
fixed through m^rn^, fn, m r)', and the scalar a . Input is indicated through 
a*. 

Table 1. 



The main parameters of the model: current quark masses m u , m s , and corresponding 
constituent masses, M u and M s in MeV, couplings G (in GeV~ 2 ) and k (in GeV -5 ), and 
two cutoffs A, A in GeV. The values of condensates are given in MeV. 





m u 


m s 


M u 


M s 


G 


— K 


A 


A 


— {uu) 


-<-> 1/3 


a 


6.3 


194 


398 


588 


13.5 


1300* 


0.82 


0* 


229 


172 


b 


6.3 


194 


398 


588 


13.5 


1300* 


0.82 


1.8* 


229 


172 


c 


6.3 


194 


398 


588 


13.4 


1370* 


0.82 


0* 


229 


172 


d 


6.3 


194 


398 


588 


11.8 


1370* 


0.82 


1.7* 


229 


172 


c 


2.8 


92 


216 


385 


3.14 


120 


1.37 


0* 


302 


314 


f 


2.1 


69 


196 


354 


2.15 


53 


1.64 


1.9* 


333 


363 



Table 2. 
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The main characteristics of the light pseudo-scalar mesons in MeV. The singlet-octet mixing 



angL 


3 Op is g 


iven in 


degrees 














m K 


u 


fK 


m v 


m v ' 


6 P 


a 


138* 


494* 


92* 


114* 


476 


986 


-14 


b 


138* 


494* 


92* 


114* 


487 


958 


-15 


c 


138* 


494* 


92* 


114* 


480 


1020 


-13 


d 


138* 


494* 


92* 


114* 


472 


959 


-15 


e 


138* 


494* 


92* 


129 


533 


1097* 


-1.2 


f 


138* 


494* 


92* 


129 


540 


1097* 


0.5 



Table 3. 

The characteristics of the light scalar nonet in MeV and the singlet-octet mixing angle S 
in degrees. 





m ao ~a (980) 


m K ^K* (800) 


m CT ~/o(600) 


77v~/o(980) 


e s 


a 


1040 


1267 


806 


1438 


24 


b 


981 


1219 


781 


1427 


24 


c 


1056 


1280 


805 


1447 


23.7 


d 


967 


1208 


762 


1426 


23.5 


e 


980* 


1029 


413 


1123 


19.5 


f 


980* 


992 


346 


1073 


18 



It is clear that the large differences observed among different pairs of sets 
come from the leading order contribution. For example the condensates, fx 
and the rj mass are strongly dependent on the value of k, which is one order of 
magnitude larger in the sets (a-d) as compared to set (e,f). This observation 
applies also to the scalar spectrum, with large changes resulting at leading 
order. They are best described in set (f), but this implies rather large values 
for fx and the condensates. 

One observes however that the corrections, although small, have the cor- 
rect trend, diminishing the splitting in the singlet-octet members of the pseu- 
doscalar and scalar spectra. Comparing sets (a,b) with (c,d), one sees that 
by enlarging the magnitude of k, the effect of the corrections get stronger, 
even for a smaller value of A. 
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We also remark the interesting fact that A cannot be arbitrarily increased, 
maintaining the remaining input on observables fixed. At values quite close 
to the ones indicated, solutions cease to exist. There is an intrinsic constraint 
on the size of corrections. 

One might be struck by the large variation in the values of k and their re- 
lation to the convergence of the loop expansion This can be understood 
by identifying the dimensionless expansion parameter of the series. Follow- 
ing 18J, where standard methods are used to justify the stationary phase 
approach to functional integrals, one obtains the dimensionless parameter 



Here the group structure factor 3/16 of $ a b c in eq.([T2|) as well as the factor 
1/3! appearing at each order in (|16|) have been taken into account. Note 
also that each term of the expansion carries a further suppression factor 1/n! 
In the present case n = 2 so that one get for sets (b), (d), (f) that ( 2 /2 is 
0.0096, 0.021, 0.0023 respectively. This attests for a fast convergence of the 
series. 

For a comparison with empirical values, we take from [27J: 



for the masses in the low lying pseudoscalar sector. The weak decay constants 
fl = 130.7 ± 0.1 ± 0.36 MeV, f e K = 159.8 ± 1.4 ± 0.44 MeV relate to ours 
through a normalization factor, thus /„- ~ 92.4 MeV and fx — 113 MeV. 
The low- lying scalar masses are presently 5 : 



5 As several data sets are presented in 27 please consult it for details. Here we indicate 
the lowest and the highest values collected from all samples. 




(55) 



m u = 1.5 -=-4 MeV, 

m d = 4 + 8 MeV, 

m s = 80 ^ 130 MeV, 

m^± = 139.57018 ± 0.00035 MeV, 

m K± = 493.677 ± 0.016 MeV, 

m v = 547 ±0.12 MeV, 

m v > = 957.78 ± 0.14 MeV, 



(56) 



mao^) =984.7 ±1.2 MeV, 
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m /o(600) = 400 + 1200 MeV, 
m Mm = 980 ± 10 MeV, 

m K *(soo) = 701 970 MeV. (57) 

The recent update of the light-quark condensate is ((uu+dd) /2)(1 GeV) = 
— (242 ± 15 MeV) 3 , the flavour breaking ratio is known to be (ss) / ((uu + 
dd)/2) = 0.8 ±0.3 |2EI. 

Finally, as our numerical calculations do not differ significantly from the 
leading order values, we refer to J2| where a thorough discussion of our 
leading order results has been done in comparison with the ones obtained 
from other approaches [2331 -|H2|. 



5. Conclusions 

The bosonization of the model combining the NJL and the 't Hooft multi- 
quark interactions leads to corrections associated with the stationary phase 
integration over auxiliary bosonic variables in the functional integral of the 
theory. The purpose of the present work has been to quantify NLO correc- 
tions, and to study their phenomenological effect on the mass spectrum of 
light pseudoscalar and scalar mesons. 

To this end the first correction to the tree level effective action has been 
considered. We have obtained the linear and quadratic terms (in the external 
mesonic fields) of the NLO Lagrangian. The group structure of the SU (2) x 
U(l) flavour symmetry considered leads to quite intricate expressions for the 
mass corrections. We have shown in a transparent way that they comply 
with the QCD low energy theorems. We have calculated the mass spectra of 
the low lying pseudoscalars and scalars, quark condensates and weak decay 
constants f n , fx- The corrections are small and improve slightly the leading 
order results. 

We conclude from these calculations that the series considered is well 
convergent. It is an important conclusion, because it justifies the leading 
order estimates made before, on one hand, and reports on the selfconsistency 
of the stationary phase approach applied to the bosonization of effective 
multi-quark interactions, on the other hand. 

At the same time one may still expect some noticeable effects of the NLO 
terms which can show themselves in three and higher order mesonic ampli- 
tudes, especially in the cases where there is a strong cancellation between 
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the tree level contributions. 
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Appendix 

The equations and algebra leading to the coefficients of the series 
for s a and p a can be found in [T7j . The explicit expressions for the case of 
one and two lower indices are also given there. In this appendix we collect 
explicit expressions for the coefficients h^ c and h^ c entering in the expansion 
of s a and needed for the evaluation of meson mass terms. Due to the trace 
structure of (JHU|) and since L contributes only in the blockdiagonal, with 
nonvanishing entries in the diagonal and elements of 0, 8 mixing , only the 
elements of the diagonal and (0, 8) mixing of L2 will contribute to the mass 
terms. For those we need: 



hn] n = f= '—r- ttt for a G 1, 2, 3. 

0aa W6G 3 fi+(l-u a ) 2 

h9~] n = %z for a G 4, 5, 6, 7. 

0aa ' ' ' 

. (1) « l+U) s + L) u ^100 
hgaa = ~ — 7=7^ T-, ^ f O r « ^ 1, 2, 3. 

16V3G 3 - c^ s ,r 
32V3G 3 - uj u ) z 

h$ 8 = 16v/ ^3 3 (1 + 2w„)[l + u; a (l - 2u u )} 

^008 = ^080 — ^800 — g^Q3^3 U u( U u ~ W s)(l ~~ U u) 
h 808 = h 880 = 16v ^ 3 3 ( 1 + U * + ~ 4 ^«) 

t& = 16v ^ 3// 3 C 1 + w « + ^)(! + 2 ^) 2 ( 59 ) 
(2) « l+u s -2u; u foraGl23 

^ = ~ i C z^s — T\ T~! ^2 for a G 4 ' 5 ' 6 ' 7 - 
16V6G 3 + w«) 

, (2) « 1 + ^ + ^ 

^8aa = — 7^ TT" w f O r a G 1, 2, 3. 

16V3G 3 A*+(l + 
19 



K 1 + u;„ — 2a>„ 



4os = 4 2 80 = ~ 24 :y /^G 3 fj, /i 2 ^" - + w„ - 3^ 2 ) 

^800 = 24^(^3^ ^2 ^ - ^'X 1 + ^«)( 2 + 3CJ «) 

4o8 = ^so = - 7^ 7 ( 3 + 2 ^« - + - - 12 ^«)) 



48 v / 6G 3 /i+/i 2 _ 

4s8 = ~ 48v ^3 ^ 2 (1 - 2u; u )(3 + u u - uj s - Quj u uj s - 6cj 2 ) 

(61) 



where 



Ui = 16^' = 1 ± ^ a - 2w 2 . (62) 
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